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1.  INTRODUCTION 


Under  assumed  provisions  of  a Threshold  Test  Ban  Treaty  (TTBT) , 
the  yield  of  a single  nuclear  device  which  could  be  detonated  in 
an  underground  test  would  be  limited  to  150  kt  or  less.  The  terms 
of  the  assumed  treaty  also  provide  for  Peaceful  Nuclear  Explosions 
(PNE)  utilizing  multiple  nuclear  detonations  having  a greater 
combined  yield.  It  has  been  suggested  that  a party  might  attempt 
to  evade  the  terms  of  an  TTBT  by  detonating  a single  clandestine 
event  having  a yield  greater  than  150  kt  in  the  environment  of  a 
high  yield  PNE.  This  report  examines  the  seismic  aspects  of  an 
evasion  scenario.  Specifically,  recorded  ground  motion  from 
single  contained  events  are  scaled  to  an  appropriate  yield  and 
superimposed  with  the  ground  motion  recorded  from  a multiple, 
surface  cratering  event  detonated  at  the  Nevada  Test  Site  (NTS). 
The  resulting  data  are  then  analyzed  and  compared  to  determine  the 
observable  effect  of  the  clandestine  event. 

This  report  makes  the  following  assumptions  concerning  a Threshold 
Test  Ban  Treaty: 

1.  The  maximum  yield  of  a single  nuclear  device  is  150  kt. 

2.  Peaceful  Nuclear  Experiments  are  permitted  involving 
simultaneous  multiple  explosions.  The  yield  of  each 
device  is  limited  to  150  kt  or  less. 

3.  On-site  scientific  instrumentation  is  permitted  for  the 
purposes  of  yield  estimation  and  seismic  monitoring. 
The  seismic  instruments  can  be  deployed  anywhere  within 
a 20  km  radius  of  the  PNE. 

The  above  assumptions  will  be  used  to  examine  the  constraints  on 
the  parameters  associated  with  the  clandestine  event  relative  to 
the  PNE. 

In  this  report  the  PNE  evasion  scenario  will  be  simulated  using 
recorded  ground  motion  from  a low  yield  (5.5  kt)  surface  row- 
cratering event  to  represent  the  PNE.  This  event,  Buggy  I,  is  the 


only  event  available  (Cassity,  et  al.,  1969)  which  approximates, 
except  for  yield,  the  conditions  associated  with  the  PNE  in  the 
evasion  scenario.  However,  by  suitably  scaling  the  yield  associ- 
ated with  the  clandestine  event  relative  to  the  yield  of  the  PNE, 
the  resulting  simulation  should  reasonably  represent  the  evasion 
scenario. 

In  Section  2,  the  constraints  on  yield,  detonation  delay  time,  and 
location  of  the  clandestine  event  relative  to  the  PNE  are  exam- 
ined. In  Section  3,  the  seismic  characteristics  of  a multiple, 
surface  cratering  event  are  analyzed  and  compared  to  the  seismic 
characteristics  of  a single  contained  event  of  equivalent  yield. 
In  Section  4,  the  superimposed  seismic  data  are  analyzed.  Section 
5 presents  a summary  and  the  conclusions  of  this  study  and  makes 
recommendations  for  further  studies. 
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2.  EXAMINATION  OF  THE  CONSTRAINTS 

ON  THE  YIELD,  LOCATION,  AND  DETONATION 
TIME  OF  THE  CLANDESTINE  EVENT 

2.1  INTRODUCTION 

In  order  to  simulate  the  PNE  evasion  scenario  described  in  Section 
1,  the  parameters  of  yield,  location  and  detonation  delay  time  of 
the  clandestine  event  (CE)  relative  to  the  PNE  must  be  determined. 
The  CE  parameters  are  constrained  by  both  on-site  and  far-field 
scientific  monitoring  equipment  whose  function  is  to  provide 
information  for  the  detection  of  an  evasion.  The  purpose  of  this 
section  is  to  examine  the  constraints  of  the  CE  yield,  location, 
and  detonation  delay  time  imposed  by  the  monitoring  instrumenta- 
tion and  to  determine  the  parameters  that  will  be  used  in  the 
simulation  described  in  Section  4.  Section  2.2  deals  with  the 
question  of  yield  while  Section  2.3  considers  the  questions  of 

location  and  detonation  delay. 

2.2  CONSTRAINTS  ON  YIELD 

The  evasion  scenario  involving  a PNE  and  a CE  can  be  stated  as 

follows:  a PNE  consisting  of  n devices  would  be  announced,  each 

device  having  a yield,  W,  for  a total  PNE  yield  of  W>p  = nW.  In 
order  to  accommodate  a CE  of  yield  Wc , the  announced  yield  of  each 
device  in  the  PNE  would  have  to  be  reduced  by  a fractional  amount, 
f.  The  result  would  be  a PNE  having  a modified  total  yield  W*p  = 
nfW  and  a CE  having  a yield  Wc  = (l-f)nW.  The  sum  of  W-p  and  Wc 

would  then  equal  the  total  yield  of  the  announced  PNE. 

In  order  to  detect  a possible  reduction  in  the  total  yield  of  the 
PNE,  specific  instrumentation  (slifer  devices)  will  be  deployed  to 
monitor  the  yield  of  each  nuclear  device  in  the  PNE.  All  the 
individual  yield  estimates  for  each  nuclear  device  will  be  summed 
and  then  compared  with  the  announced  total  yield  of  the  PNE. 
Since  there  is  uncertainty  associated  with  estimating  yield,  the 
fractional  reduction  in  yield  must  be  chosen  to  be  compatible  with 
the  observational  uncertainty  in  yield  determination.  Based  on 
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present  yield  estimation  capability,  it  appears  that  a reduction 
in  yield  of  approximately  20  percent  is  reasonable  so  that  an  f = 
0.8  would  be  appropriate. 

In  addition  to  on-site  instrumentation  for  estimates  of  yield,  the 
PNE  would  also  be  monitored  by  teleseismic  instrumentation  which 
also  permits  yield  estimation.  This  yield  determination  would  be 
for  the  combined  yield  of  both  the  PNE  and  the  clandestine  device. 
Again,  the  yield  determination  will  contain  some  uncertainty. 
This  fact  would  be  used  to  increase  the  yield  of  the  clandestine 
device,  and  still  operate  within  the  yield  determination  uncer- 
tainty of  the  teleseismic  information.  If  a 20  percent  uncer- 
tainty is  assumed,  then  multiplying  the  total  yield  by  1.2  would 


give  a reasonable  value  for  the  combined  yield  of  the  PNE  and  CE. 
If  the  increased  yield  is  associated  with  the  CE,  it  is  easy  to 
show  that  the  CE  yield,  Wc,  would  be  0.4  Wp,  where  Wp  is  the 
announced  total  yield  of  the  PNE. 

Thus,  based  on  on-site  and  far-field  yield  estimation  capability, 
the  range  on  yield  for  the  CE  would  be  from  0.2  Wp  to  0.4  Wp  where 
W-p  is  the  announced  total  yield  of  the  PNE.  The  total  yield  of 
the  modified  PNE  would  be  0.8  W-p.  These  parameters  will  be  used 
in  Section  4 to  determine  the  CE  yield  relative  to  the  yield 
associated  with  the  PNE. 

2.3  CONSTRAINTS  ON  THE  CLANDESTINE  EVENT  LOCATION  AND  DETONATION 
DELAY  TIME 

The  terms  of  our  assumed  threshold  test  ban  treaty  also  provide 
for  the  deployment  of  on-site  seismic  instruments  to  monitor 
ground  motion.  The  amplitude  and  time  of  first  arrival  of  seismic 
motion  from  both  the  PNE  and  CE  are  a direct  function  of  the  rela- 
tive location  and  time  of  detonation  associated  with  the  two 
sources  of  energy.  This  section  examines  the  constraints  imposed 
on  location  and  detonation  delay  time  by  both  on-site  and  far- 
field  seismic  instrumentation. 


Each  on-site  seismograph  will  record  the  sum  or  superposition  of 
the  ground  motion  caused  by  the  PNE  and  the  CE.  For  a particular 
seismic  recording  station,  the  peak  amplitude  levels  of  the  PNE 
and  CE  wave  trains  are  controlled  by  both  the  yields  and  recording 
distances  of  the  two  sources.  The  time  associated  with  the  first 
motion  arrival  from  each  source  is  controlled  by  the  respective 
source-to-station  distances,  the  wave  propagation  velocity  in  the 
area  and  by  any  time  delay  between  the  two  detonations.  Formulas 
for  estimating  the  peak  amplitude  and  first  arrival  time  will  be 
given  below  and  will  be  then  used  to  examine  the  constraints  on 
the  location  of  the  CE  relative  to  the  PNE  based  on  these  general 
parameters . 

For  a single  contained  event  of  yield,  Wc,  and  source-to-station 
distance,  Rc , the  expected  vertical  peak  velocity,  vc,  is  given  by 
(Environmental  Research  Corp.,  1974) 

v c = 4.89  Wc°-74  rc-1.47  (2-1) 

An  empirical  formula  for  the  peak  velocity  associated  with  a multi- 
ple surface  cratering  event  is  not  available.  However  it  will  be 
shown  in  Section  3,  that  the  peak  velocity  from  an  observed  multi- 
ple detonation  is  about  a factor  of  1.5  less  than  the  peak  veloc- 
ity from  a single  contained  event.  Using  this  empirical  factor, 
the  expected  peak  velocity,  vp,  of  a PNE  of  combined  yield,  Wp, 
and  source-to-station  distance,  Rp,  is  given  by 

vp  = 3.25  Wp0.74  Rp-1.47  (2-2) 

Also,  an  empirical  equation  is  not  available  for  estimating  the 
peak  velocity  associated  with  the  superposition  of  two  seismic 
wave  trains.  However,  a reasonable  estimate  can  be  obtained  using 
a root-mean-squa re  estimate  given  by 

vpc  = Vvp2  + vc2  (2-3) 
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for  a PNE  with  a 


The  equation  for  the  time  of  first  arrival,  tp, 
source-  to-station  distance,  Rp,  is  given  by 


t 


P 


(2-4) 


when  V is  the  propagation  velocity  associated  with  P-waves.  For 
the  CE  which  is  detonated  with  a time  delay  of  Atc  seconds,  the 
time  of  first  arrival,  tc,  for  a source-to-stat ion  distance,  Rc, 
is  given  by 


R 

+ A t n 

v c 


(2-5) 


In  actual  practice,  the  propagation  velocity,  V,  would  not  be  the 
same  for  all  transmission  paths,  but  for  initial  estimation  pur- 
poses, a single  value  will  be  assumed. 


The  equations  developed  above  will  now  be  applied  to  the  case  in 
which  the  CE  is  located  outside  the  on-site  seismic  array.  Assume 
an  announced  PNE  of  1500  kt,  a modified  PNE  of  1200  kt  and  a CE  of 
600  kt  located  40  km  from  the  PNE.  Consider  two  seismic  sta- 
tions on  a line  connecting  the  two  sources,  located  at  10  and  20 
km  from  the  PNE.  Assume  also  that  the  body  wave  propagation 
velocity  is  4 km/sec  and  that  the  events  are  detonated  simultane- 
ously. For  the  station  located  at  10  km,  the  expected  peak  veloci- 
ties and  first  arrival  times  for  the  individual  and  superimposed 
wave  trains  are  as  follows: 


Announced  PNE: 
Modified  PNE: 

CE: 

Modified  PNE+CE: 


Expected 
Peak  Velocity 

24.7  cm/sec 
20.9  cm/sec 
3.75  cm/sec 
21.2  cm/sec 


First 

Arrival  Time 

2.5  sec 

2.5  sec 

7.5  sec 

2.5  sec 


These  data  show  that  the  CE  peak  velocity  would  be  a factor  of 
about  5.6  below  that  of  the  PNE.  The  time  of  first  arrival  would 
appear  normal.  The  expected  peak  velocity  from  the  announced  PNE 
and  the  modified  PNE  + CE  differ  by  approximately  17  percent,  a 
difference  which  is  small  compared  to  normal  variation  of  observed 
peak  velocity  about  the  expected  mean  value.  Thus,  for  the  sta- 
tion at  10  km,  the  peak  velocity  and  first  arrival  time  from  the 
modified  PNE  + CE  would  differ  little  from  that  expected  for  the 
announced  PNE. 


For  the  station  located  20  km  from  both  the  PNE  and  the  CE,  the 
expected  peak  velocities  and  first  arrival  times  are  as  follows: 


Announced  PNE 
Modified  PNE: 

CE: 

Modified  PNE+CE: 


Expected 
Peak  Velocity 

8.90  cm/sec 
7.55  cm/sec 
6.80  cm/sec 
10.16  cm./sec 


First 

Arrival  Time 

5 sec 
5 sec 
5 sec 
5 sec 


In  this  case  the  peak  velocities  from  the  modified  PNE  and  CE  are 
nearly  equal.  The  expected  peak  velocity  from  the  modified  PNE+CE 
is  approximately  14  percent  greater  than  that  expected  from  the 
announced  PNE,  again  an  insignificant  variation.  In  all  cases  the 
time  of  first  arrival  would  be  the  same. 


For  the  two  cases  considered  above,  it  appears  that  based  on  the 
results  of  near-field  instrumentation,  the  modified  PNE+CE  would 
be  indistinguishable  from  the  announced  PNE  and  that  a source 
separation  of  40  km  would  be  a viable  evasion  tactic.  However, 
analysis  of  teleseismic  data  would  lead  to  the  conclusion  that  a 
CE  had  been  detonated,  as  the  following  discussion  will  show. 

Considerable  effort  has  been  expended  in  a program  to  identify  two 
distinct  sources  of  seismic  energy  based  on  teleseismic  data.  A 
primary  factor  in  distinguishing  two  sources  is  the  time  separa- 
ting T,  between  the  initial  arrivals.  The  time  separation  is  a 
function  of  both  the  distance,  AD,  between  the  sources,  the  detona- 
tion delay  time.  At,  and  the  distance.  A,  to  the  seismic  recording 
station.  The  present  discrimination  capability  is  shown  in  Figure 
2.1  (Alewine,  1978),  which  expresses  the  probability  of  making  a 
specified  error  in  classification  as  a function  of  source  sepa- 
ration and  detonation  delay  for  recording  stations  at  distances 
A=30°  and  60°. 


The  information  expressed  in  Figure  2.1  can  be  stated  in  a mathe- 
matical equation.  For  example,  equations  for  A = 60°  and  for 

class i f icat ion  errors  of  15  percent  and  50  percent  are  given  by 

P(E)  = 0.15;  0.0625  AD  + At  = 2.5  sec  (2-6) 

P(E)  = 0.50;  0.0625  AD  + At  = 0.6  sec  (2-7) 

These  equations  can  be  considered  as  constraints  on  AD  and  At  for 
the  specified  teleseismic  detection  capability. 


* 
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FIGURE  2.1:  TELESEISMIC  DISCRIMINATION  OF  TWO  SEISMIC  SOURCES  BASED  ON  SOURCE 
SEPARATION,  DETONATION  DELAY  TIME  AND  TELESEISMIC  DISTANCE 


Returning  now  to  the  PNE  and  CE  evasion  scenario  in  which  the 
location  separation  is  40  km  with  zero  time  delay,  the  previous 
analysis,  based  on  near-field  seismic  data,  indicated  that  it 
would  be  difficult  to  identify  the  CE  based  on  peak  velocity  and 
first  arrival  time.  However,  this  case  would  be  detectable  with  a 
15  percent  probability  of  error  based  on  teleseismic  information. 
In  fact,  the  greater  the  PNE  and  CE  source  separation,  the  smaller 
the  error  in  discriminating  the  two  sources  from  teleseismic  data. 
At  the  same  time,  the  greater  the  difficulty  in  detecting  the 
clandestine  event  based  on  on-site  seismic  information. 

Examination  of  Figure  2.1  shows  that  as  the  source  separation 
decreases,  the  error  in  discrimination  capability  increases  for  a 
given  detonation  delay.  But  as  the  source  separation  decreases, 
the  clandestine  event  comes  closer  to  the  on-site  seismic  array. 
The  consequences  in  terms  of  peak  velocity  and  first  arrival  time 
is  shown  in  the  following  example. 

Assume  that  the  source  separation  is  30  km  and  that  a seismic 
station  is  located  20  km  from  the  PNE  and  10  km  from  the  CE.  Let 
the  teleseismic  discrimination  error  be  15  percent,  so  that  the 
permissible  detonation  delay,  is  At  = 0.625  seconds.  The  expected 
peak  velocities  and  first  arrival  times  at  the  seismic  station  are 
as  follows: 


Announced  PNE 
Modified  PNE 
CE 

Modified  PNE+CE 


Expected 
Peak  Velocity 

8.90  cm/sec 
7.55  cm/sec 
18.84  cm/sec 
20.29  cm/sec 


First 

Arrival  Time 

5.0  sec 
5.0  sec 
3.12  sec 
3.12  sec 


These  data  show  that  the  expected  velocity  from  the  modified 
PNE+CE  is  more  than  a factor  of  two  greater  than  would  be  expected 
from  the  announced  PNE.  In  addition,  the  first  arrival  time  would 
be  1.88  seconds  sooner  than  expected.  These  data  would  certainly 
appear  anomalous,  especially  the  significantly  smaller  first 
arrival  time. 


As  the  source  separation  becomes  less  than  30  km,  the  situation 
becomes  worse  from  the  on-site  seismic  viewpoint.  It  seems  cer- 
tain that  the  clandestine  event  would  not  be  located  near  possible 
seismic  station  locations.  The  one  location  for  the  CE  which  does 
appear  feasible  is  directly  below  the  PNE,  i.e.,  zero  surface  sepa- 
ration. This  location  would  also  be  desirable  from  a teleseismic 
discrimination  viewpoint.  If  the  detonation  delay  were  restricted 
to  0.3  seconds  or  less,  Figure  2.1  indicates  that  the  probability 
of  error  in  discrimination  would  be  50  percent  or  greater  based  on 
teleseismic  data. 


For  zero  separation  and  a 0.3  second  detonation  delay,  the 
expected  peak  velocities  and  first  arrival  times  for  an  on-site 
seismic  station  located  at  10  km  are: 


Announced  PNE 
Modified  PNE 
CE 

Modified  PNE+CE 


Expected 
Peak  Velocity 

24.7  cm/sec 
20.9  cm/sec 
18.89  cm/sec 
28.1  cm/sec 


First 

Arrival  Time 

2 . 5  sec 

2.5  sec 
2.8  sec 

2.5  sec 


And  for  a station  located  at  20  km  they  are: 


Announced  PNE 
Modified  PNE 
CE 

Modified  PNE+CE 


Expected 
Peak  Velocity 

8.90  cm/sec 
7.55  cm/sec 
6.80  cm/sec 
10.16  cm/sec 


First 

Arrival  Time 

5.0  sec 

5.0  sec 
5 . 3 sec 

5.0  sec 


The  above  data  show  that  the  expected  peak  velocity  from  the  modi- 
fied PNE+CE  would  be  only  slightly  greater  than  expected  from  the 
announced  PNE.  In  addition,  the  time  of  first  arrival  would 
appear  normal. 


Based  on  the  above  analysis,  it  is  concluded  that  the  most  likely 
location  for  the  clandestine  event  would  be  directly  below  the 
PNE.  A detonation  delay  of  0.3  seconds  or  less  would  minimize  the 
possibility  of  distinguishing  the  two  sources  of  energy  based  on 
an  analysis  of  teleseismic  data. 


3.  COMPARISON  OF  A SURFACE  ROW-CRATERING  EVENT  WITH 
SINGLE  CONTAINED  EVENTS  OF  EQUIVALENT  YIELD 

3.1  INTRODUCTION 

In  simulating  the  ground  motion  for  a PNE  evasion  scenario, 
recorded  ground  motion  from  a single  contained  underground  nuclear 
explosion  will  be  superimposed  with  recorded  ground  motion  from  a 
multiple,  surface  row-cratering  event  representing  the  PNE. 
Because  of  the  marked  difference  in  the  source  configuration  of 
these  two  classes  of  events,  the  characteristics  of  the  seismic 
waves  from  the  two  types  of  sources  should  also  be  different. 
This  section  examines  quantitative  differences  in  peak  amplitude 
and  spectral  characteristics  of  ground  motion  from  these  two 
classes  of  events.  In  addition,  this  section  describes  the 
available  seismic  data  sample  that  will  be  used  in  the  simulated 
PNE  evasion  scenario. 

3.2  THE  SURFACE  ROW-CRATERING  EVENT 

On  March  12,  1968,  a plowshare  experiment,  Buggy  I,  was  conducted 
at  the  Nevada  Test  Site  and  involved  the  near  surface  detonation 
of  five  nuclear  devices  with  a combined  yield  of  5.5  kt  (Cassity, 
et  al.,  1969).  The  nuclear  devices  had  yields  of  1.1  kt  and  were 
buried  at  a depth  of  135  ft.  The  devices  were  in-line  with  a hori- 
zontal spacing  of  135  ft.  which  gives  a 600-foot  spacing  between 
the  end  devices.  Buggy  I is  the  only  surface  row-cratering  event 
for  which  seismic  recordings  are  available. 

The  available  seismic  data  from  the  Buggy  I event  consists  of  the 
vertical  component  of  ground  velocity  recorded  at  four  stations. 
The  stations  are  designated  as  447,  458,  449,  and  450  and  are 
located,  respectively,  at  surface  distances  of  2.69,  2.98,  6.37, 
and  11. 1 km  from  the  center  of  the  Buggy  I configuration.  The 
relationship  of  each  seismic  station  to  Buggy  I is  given  in  Figure 
3.1  which  shows  three  stations  in-line  and  one  station  at  135° 
relative  to  the  direction  defined  by  the  nuclear  devices. 
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FIGURE  3.1:  THE  LOCATION  OF  THE  SEISMIC  RECORDING  STATIONS 
RELATIVE  TO  THE  BUGGY  I EVENT 
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The  observed  motion  from  Buggy  I will  be  used  as  the  PNE  in  the 
evasion  scenario.  In  Section  2,  it  was  shown  that  the  yield  of 
the  CE  could  range  between  20%  to  40%  of  the  announced  PNE  yield 
and  still  fall  within  the  yield  determination  uncertainty  of  both 
on-site  and  far-field  instrumental ion.  For  the  Buggy  I total 
yield  of  5.5  kt,  which  would  be  the  modified  PNE  yield  or  80%  of 
the  announced  PNE  yield,  the  yield  association  with  the  announced 
PNE  would  be  5.5  kt/0.8  = 6.875  kt.  The  yield  range  association 
with  the  CE  would  be  20%  to  40%  of  the  announced  PNE  yield  or  a 
yield  range  of  1.375  kt  to  2.75  kt.  In  Section  4,  simulations 
will  be  performed  for  two  cases  for  the  CE  yield.  The  actual 
yield  values  used  are  1.375  kt  and  2.5  kt. 


3.3  SINGLE  CONTAINED  EVENTS 

A search  of  the  available  data  bank  of  NTS  seismograms  revealed  a 
sample  of  26  vertical  velocity  recordings  at  stations  within  12  km 
of  the  detonation  point.  Within  this  set,  11  matched  the  source- 
to-station  distances  associated  with  the  Buggy  I event.  None  of 
the  events  associated  with  the  velocity  seismograms  had  yields  of 
1.375  or  2.5  kt  that  are  required  for  the  PNE  evasion  simulation. 

Fortunately,  procedures  have  been  developed  to  scale  observed 
seismogram  time  histories  for  the  source  parameters  of  yield  and 
depth  of  burial  (Mueller  and  Murphy,  1971;  Murphy,  1977).  Essen- 
tially, an  observed  time  history  is  operated  on  by  a frequency- 
dependent  transfer  function  containing,  as  parameters,  the  actual 
and  desired  yield  and  depth  of  burial.  The  result  is  an  estimated 
tine  history  scaled  to  the  desired  yield  and  depth  of  burial. 
Examples  of  transformed  seismograms  are  given  in  Murphy,  1977. 
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3.4  COMPARISON  OF  THE  SEISMIC  CHARACTERISTICS  OF  A SURFACE  ROW- 
CRATERING AND  SINGLE  CONTAINED  EVENTS  OF  EQUIVALENT  YIELD 

The  seismic  characteristics  of  time  histories  from  single  con- 
tained events  detonated  at  the  NTS  have  been  extensively  studied 
since  the  beginning  of  the  underground  nuclear  testing  program. 

Based  on  analyses  of  observed  data,  empirical  equations  have  been 
developed  to  estimate  the  peak  amplitude  and  spectral  amplitude  of 
ground  motion  as  a function  of  the  yield  of  the  device  and  the 
source- to-stat ion  distance.  Because  of  the  limited  experience 
with  surface  row-cratering  events,  amplitude  and  spectral  scaling 
relationships  are  unknown.  The  purpose  of  this  section  is  to 
compare  the  peak  amplitude  and  spectral  characteristics  of  Buggy  I 
seismograms  with  those  characteristics  obtained  from  a set  of 
close-in  seismograms  from  single  contained  events  scaled  to  the 
equivalent  yield. 

( 

The  peak  velocities  from  single  contained  events  scaled  to  5.5  kt 
have  been  analyzed  statistically  using  a power  law  relationship 
given  by 

vp  = v0R-n  (3-1) 

where  Vp  and  R are  observations  and  vg  and  n are  constants  to  be 
determined.  An  analysis  results  in  the  equation 

Vp  = 10. 9R-1*67  (3-2) 

where  Vp  is  the  expected  peak  velocity  in  cm/sec  and  the  source- 
to-station  distance,  R,  is  in  km.  The  statistical  analysis  also 
gives  an  estimate  of  the  scatter  of  the  observations  about  the 
expected  value.  This  measure  of  scatter  is  called  the  standard 
error  of  estimate,  o , and  the  analysis  gives  a = 2.1.  Multiplying 
and  dividing  the  expected  peak  velocity,  Vp,  by  defines  an 
interval  in  which  67%  of  the  observations  are  expected  to  fall. 


A plot  of  equation  3.1  together  with  the  observed  peak  velocities 
from  the  scaled  single  contained  events  and  the  observed  peak 
velocities  from  Buggy  I are  shown  in  Figure  3.2.  As  this  figure 
indicates,  the  peak  velocities  from  Buggy  I fall  below  the  line 
defining  the  expected  peak  velocity.  On  the  average,  the  Buggy  I 
peaks  are  a factor  of  1.5  below  the  expected  value  for  a single 
contained  event  of  equivalent  yield.  Figure  3.2  also  shows  the 
considerable  scatter  associated  with  observed  peak  velocities.  In 
addition  to  peak  velocity,  the  single  contained  events  have  been 
analyzed  to  determine  spectral  information  using  a narrow-band 
filter  technique.  The  filters  utilized  in  the  analysis  are  second 
order  ac  ive  filters  characterized  by  a center  frequency,  f^,  and 

a damping  factor,  h,  which  was  chosen  to  be  5 percent  of  critical 

damping.  The  transfer  function,  F(S),  of  this  filter  is  given  by 

2hw.  2 

F(S)  = — ± 2 ( 3-3  ) 

S ^ + 2hw.S  + w. 

l l 

where  S is  the  Laplace  transform  variable,  and  w^  = 2nf^.  Each 
seismogram  is  processed  through  a set  of  narrow-band  filters  each 
having  a different  center  frequency.  The  peak  amplitude  response 
from  each  filter  represents  the  spectral  value  Svi  for  each  center 
frequency,  f^.  it  is  customary  to  plot,  Svi,  versus  period,  p^, 
which  is  related  to  frequency  by  p^  = 1/fi.  An  example  of  the 

spectral  information  obtained  using  this  technique  is  shown  in 

Figure  3.3. 

As  was  done  with  peak  velocity,  spectral  information  was  analyzed 
statistically  using  a power  law  given  by 

svi  = SvoiR~ni  i = 1,  2...N.  (3-4) 

where  Sv^  and  R are  observed  and  Svoi  and  n^  are  parameters  to  be 
estimated  for  each  period,  pj.  The  analysis  results  which  give 
estimates  of  Svo£,  nj  and  the  standard  error  of  estimate,  sj,  are 
shown  in  Table  3.1 
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nen.! RE  3.2:  A PLOT  OF  THE  PEAK  VELOCITY  VERSUS  DISTANCE 
FOR  BUGGY  I (CIRCLES)  AND  SINGLE  CONTAINED 
EVENTS  (DOTS) 


A plot  of  the  expected  spectral  value  for  the  period  0.499  seconds 
is  shown  in  Figure  3.4  together  with  the  observed  spectral  values 
from  single  contained  events  and  from  the  Buggy  I event.  These 
data  show  that  the  spectral  values  from  Buggy  I are  all  lower  than 
that  expected  from  a single  contained  event.  The  average  ratio  of 
expected  to  observed  Buggy  I is  2.5  at  this  period.  The  average 
ratios  for  the  13  periods  examined  are  shown  in  Figure  3.5.  At 
all  but  one  period,  the  ratio  is  greater  than  one. 

In  terms  of  both  peak  amplitude  and  spectral  values,  single  con- 
tained events  can  be  expected  to  give  greater  amplitudes  than  a 
surface  row-cratering  event  of  equivalent  yield. 

Section  3 has  examined  the  differences  in  peak  amplitude  and  the 
spectral  characteristics  of  the  Buggy  I surface  row-cratering 
event  and  single  contained  events  of  equivalent  yield.  In  Section 
4,  ground  motion  from  these  two  classes  of  events  will  be  superim- 
posed to  simulate  the  PNE  evasion  scenario  and  will  be  analyzed  to 
determine  the  effects  caused  by  the  clandestine  event. 
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ESTIMATED  PARAME 
OF  SPECTRAL  INFORMATION 


Period 

S 

VO 

0.050 

2 .95 

0.067 

3.98 

0.100 

9.54 

0.150 

26.0 

0 .224 

52.8 

0.334 

29.1 

0.499 

25.2 

0.743 

18.67 

1.110 

12.83 

1.650 

3.90 

2.470 

1 .10 

3.680 

1 . 368 

6 . 070 

. 182 

m 

a 

s 

2.03 

2.19 

1.95 

2.22 

2.13 

2.21 

2.48 

2.13 

2.58 

2.33 

1.72 

2.21 

1.58 

2.13 

1 . 58 

2.49 

1.48 

2.72 

0.98 

2.46 

0 .72 

2.18 

1 . 26 

2.71 

0 .135 

2.35 

4.  SIMULATION  AND  ANALYSIS  OF 
THE  PNE  EVASION  SCENARIO 

' 

4 . 1 INTRODUCTION 

The  purpose  of  this  section  is  to  present  the  results  of  a simu- 
lated PNE  evasion  scenario  using  observed  seismograms  from  the 
Buggy  I row-cratering  events  and  single  contained  events  scaled  to 
two  simulation  yields:  1.375  kt  and  2.5  kt.  Seismograms  from 

these  two  classes  of  events  are  summed  using  a time  delay  of  0.3 
seconds.  Spectral  content  of  the  Buggy  I seismograms  will  be 
compared  to  the  spectral  content  of  the  superposition. 

i 

4.2  DATA  SAMPLE 

Seismograms  from  the  single  contained  event  data  sample  were 
selected  on  the  basis  of  the  source-to-station  distances  associ- 
ated with  Buggy  I.  General  information  concerning  the  selected 
events  is  given  in  Table  4.1  which  identifies  each  event  and 
station  associated  with  the  Buggy  I stations  and  provides  informa- 
tion on  the  yield  and  distance  for  the  single  continued  events. 
Seven  of  the  13  seismograms  are  associated  with  the  Almendro  event 
which  has  the  highest  yield  in  the  sample.  Almendro  stations  L03 
and  L04  are  used  to  form  superpositions  with  Buggy  I stations  447 
and  458,  even  though  the  source-to-station  distances  are  only 
approximately  the  same. 

Each  of  the  single  contained  event  seismograms  was  scaled  to  the 
CE  yields  of  1.375  and  2.5  kt  using  the  source  scaling  technique 
referenced  in  Section  3.  The  data  were  checked  for  reasonableness 
by  comparing  the  peak  velocity  of  the  scaled  seismograms  with  the 
peak  velocity  associated  with  Buggy  I.  Examination  of  these  data 
indicated  unreasonably  large  peak  velocities  for  some  Almendro  sta- 
tions. This  can  be  seen  in  Table  4.2  which  gives  the  peak  veloci- 
ties of  all  the  Buggy  I stations  and  the  1.375  kt  scaled  Almendro 

% 


TABLE  4.1 

GENERAL  EVENT  INFORMATION 


Buggy  I(5.5kt) 

Station 

Distance 

44  7 

2.69 

i 

458 

2.97 

449 

6.  37 

450 

11.1 

stations.  Almendro  stations  L02  and  A19T  are  respectively  a fac- 
tor of  3.5  and  4.25  greater  than  the  observed  Buggy  I peak  veloci- 
ty at  station  449.  Because  of  the  large  yield  extrapolation 
involved  in  the  source  scaling,  >200  to  1.375  kt,  these  data  may 
be  questionable.  Therefore,  an  amplitude  adjustment  factor  was 
applied  to  the  scaled  Almendro  seismograms  so  that  the  ratio  of 
observer  Buggy  I to  scaled  Almendro  agreed  with  what  would  be 
prediction  based  on  the  peak  amplitude  prediction  equations. 

4.3  SUPERPOSITION  OF  SCALED  ALMENDRO  AND  BUGGY  I 

Using  the  0.3-second  time  delay,  the  scaled  single  contained 
Almendro  event  seismograms  were  added  to  the  appropriate  Buggy  I 
seismograms  to  form  the  superposition  seismograms.  Examples 
showing  each  individual  seismogram  and  the  superposition  are  given 
in  Figures  4.1  through  4.8.  The  amplitude  scale  is  the  same  for 
each  seismogram  in  a given  figure.  Examination  of  these  figures 
indicates  that  the  scaled  Almendro  seismograms  are  relatively 
short  transient  waveforms  in  comparison  to  the  Buggy  I seismo- 
grams; consequently,  the  superposition  seismograms  differ  primar- 
ily in  the  front  part  of  the  record. 

Spectra  from  both  Buggy  I and  the  superposition  seismograms  were 
generated  and  are  shown  in  Figures  4.9  through  4.16  for  clandes- 
tine event  yields  of  1.375  and  2.5  kt.  These  data  show  that,  for 
periods  greater  than  0.6  seconds,  the  Buggy  I and  the  superposi- 
tion spectra  are  very  nearly  equal  at  all  stations.  For  periods 
shorter  than  0.6  seconds,  the  superposition  spectra  are,  in  gen- 
eral, greater  than  the  Buggy  I spectra. 

Spectral  ratios  were  computed  and  then  geometrically  averaged  to 
define  quantitatively  the  difference  between  Buggy  I and  the  super- 
positions. The  result  are  shown  in  Figure  4.17  for  the  two  yields 
used  in  scaling  the  Almendro  data.  The  maximum  ratio  occurs  at 
0.1  seconds  (10  Hz). 


FIGURE  4.2:  SEISMOGRAMS  FROM  BUGGY  I,  STATION  447  (Top  Trace);  SCALED  ALMENDRO 
(1.375  kt) , STATION  L04  (Middle  Trace);  AND  THE  SUPERPOSITION 
(Bottom  Trace) 


FIGURE  4.3:  SEISMOGRAMS  FROM  BUGGY  I,  STATION  458  (Top  Trace);  SCALED  ALMENDRO 
(1.375  kt) , STATION  L03  (Middle  Trace);  AND  THE  SUPERPOSITION 
(Bottom  Trace) 


FIGURE  4.4:  SEISMOGRAMS  FROM  BUGGY  I,  STATION  458  (Top  Trace);  SCALED  ALMENDRO 
(1.375  kt) , STATION  L04  (Middle  Trace);  AND  SUPERPOSITION 
(Bottom  Trace) 


FIGURE  4.5:  SEISMOGRAMS  FROM  BUGGY  I,  STATION  449  (Top  Trace ); SCALED  ALMENDRO 
(1.375  kt) ; STATION  L05  (Middle  Trace);  AND  THE  SUPERPOSITION 
(Bottom  Trace) 
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FIGURE  4.7:  SEISMOGRAMS  FROM  BUGGY  I.  STATION  449  (Top  Trace);  SCALED  ALMENDRO 
(1.375  kt);  STATION  A19T  (Middle  Trace);  AND  THE  SUPERPOSITION 
(Bottom  Trace) 


(1.375  kt)  STATION  L06  (Middle  Trace);  AND  THE  SUPERPOSITION 
(Bottom  Trace) 
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FIGURE 


.11:  SPECTRA  FROM  BUGGY  I,  STATION  458  (DASHED 

CURVE)  AND  THE  SUPERPOSITION  WITH  SCALED 
ALMENDRO,  STATION  L03,  1.375  kt  (SOLID 
CURVE)  AND  2.5  kt  (DOTTED  CURVE) 
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FIGURE  4.12:  SPECTRA  FROM  BUGGY  I,  STATION  458  (DASHED 
CURVE)  AND  THE  SUPERPOSITION  WITH  SCALED 
ALMENDRO,  STATION  L04 , 1.375  kt  (SOLID 
CURVE)  AND  2.5  kt  (DOTTED  CURVE) 
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FIGURE  4.1 


0 1 
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5:  SPECTRA  FROM  BUGGY  I,  STATION  449  (DASHED 

CURVE)  AND  THE  SUPERPOSITION  WITH  SCALED 
ALMENDRO,  STATION  A19T,  1.375  kt  (SOLID 
CURVE)  AND  2.5  kt  (DOTTED  CURVE) 
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4.4  SUPERPOSITIONS  INVOLVING  SEISMOGRAMS  FROM  FOUR  EVENTS 


In  addition  to  the  Almendro  Event,  four  NTS  events  (see  Table  4.1) 
had  at  least  one  seismic  recording  station  which  approximated  the 
source-to-station  distance  associated  with  a Buggy  I station.  The 
yields  of  these  events  ranged  from  less  than  one  kt  to  slightly 
more  than  50  kt,  thus  requiring  less  yield  scaling  than  Almendro 
(>200  kt)  to  obtain  the  1.375  and  2.5  kt  yields  necessary  to  simu- 
late the  evasion  scenario. 

The  individual  and  superposition  seismograms  associated  with  this 
data  sample  are  shown  in  Figures  4.18  through  4.23  for  a CE  yield 
of  2.5  kt.  In  general,  the  duration  of  the  single  contained  event 
seismograms  are  comparable  to  or  have  greater  duration  than  the 
Buggy  I seismograms,  in  contrast  to  the  relatively  short  duration 
associated  with  the  scaled  Almendro  seismograms. 

Spectral  data  from  the  Buggy  I and  superposition  seismograms  are 
shown  in  Figures  4.24  through  4.29  for  the  1.375  and  2.5  kt  yields 
of  the  clandestine  event.  Figures  4.24,  4.26,  4.27  and  4.29  show 
considerable  differences  between  the  superposition  and  Buggy  I 
spectra  while  figures  4.25  and  4.28  show  only  slight  differences. 

Spectral  ratio  information  for  this  data  set  is  shown  in  Figure 
4.30.  It  can  be  seen  that  the  spectral  content  of  the  superposi- 
tion is  greater,  on  the  average,  over  the  entire  period  range 
examined.  The  spectral  ratio  is  slightly  greater  than  a factor  of 
two  at  periods  of  0.15  seconds  (6.7  Hz)  and  1.1  seconds  (0.9a  Hz) 
for  a clandestine  yield  of  2.5  kt. 

Examination  of  the  spectral  ratio  from  both  the  Almendro  and  the 
four  additional  events  used  in  the  superpositions  indicated  an 
increase  in  spectral  amplitude  relative  to  the  spectral  amplitude 
associated  with  Buggy  I. 
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FIGURE  4.18:  SEISMOGRAMS  FROM  BUGGY  I,  STATION  447  (Top  Trace);  SCALED  EVENT  A 

(2.5  kt) , STATION  790  (Middle  Trace);  AND  THE  SUPERPOSITION  (Bottom 


FIGURE  4.19:  SEISMOGRAMS  FROM  BUGGY  I,  STATION  458  (Top  Trace);  SCALED  EVENT 
(2.5  kt) , STATION  750  (Middle  Trace);  AND  THE  SUPERPOSITION 
(Bottom  Trace) 


FIGURE  4.20:  SEISMOGRAMS  FROM  BUGGY  I,  STATION  458  (Top  Trace);  SCALED  EVENT  A 

(2.5  kt) , STATION  780  (Middle  Trace);  AND  THE  SUPERPOSITION  (Bottom 


FIGURE  4.21:  SEISMOGRAMS  FROM  BUGGY  I,  STATION  458  (Top  Trace);  SCALED  EVENT 
(2.5  kt) , STATION  L08  (Middle  Trace);  AND  THE  SUPERPOSITION 
(Bottom  Trace) 


FIGURE  4.22:  SEISMOGRAMS  FROM  BUGGY  I,  STATION  449  (Top  Trace);  SCALED  EVENT  C 

(2.5  kt) , STATION  765  (Middle  Trace);  AND  THE  SUPERPOSITION  (Bottom 
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FIGURE 
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.24:  SPECTRA  FROM  BUGGY  I,  STATION  447  (DASHED 

CURVE) , AND  THE  SUPERPOSITION  WITH  SCALED 
EVENT  A,  STATION  790,  1.375  kt  (SOLID 
CURVE)  AND  2.5  kt  (DOTTED  CURVE) 
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FIGURE 


.27:  SPECTRA  FROM  BUGGY  I,  STATION  458  (DASHED 

CURVE)  AND  THE  SUPERPOSITION  WITH  SCALED 
EVENT  B,  STATION  L08 , 1.375  kt  (SOLID  CURVE) 
AND  2.5  kt  (DOTTED  CURVE) 
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FIGURE 


.28:  SPECTRA  FROM  BUGGY  I,  STATION  449  (DASHED 

CURVE)  AND  THE  SUPERPOSITION  WITH  SCALED 
EVENT  C,  STATION  765,  1.375  kt  (SOLID  CURVE) 
AND  2.5  kt  (DOTTED  CURVE) 
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5.  SUMMARY,  CONCLUSIONS  AND  RECOMMENDATIONS 


5.1  SUMMARY  AND  CONCLUSIONS 

This  report  has  examined  the  seismic  aspects  of  a simulated  PNE 
evasion  scenario  in  which  a single  contained  nuclear  device  is 
detonated  in  a multiple  surface  row-cratering  environment.  The 
constraints  on  the  yield,  detonation  delay  time  and  location  of 
the  clandestine  event  relative  to  the  PNE  are  examined  in  terms  of 
the  detection  capability  of  both  far-field  and  on-site  data.  This 
examination  indicated  that  the  yield  of  the  clandestine  event 
could  range  from  20%  to  40%  of  the  announced  PNE  yield.  In  addi- 
tion, the  examination  indicated  that  the  most  probable  location  of 
the  clandestine  event  would  be  directed  below  the  PNE  location  and 
that  the  detonation  delay  would  range  from  zero  to  0.3  seconds. 
The  spectral  characteristics  of  seismic  data  from  a multiple 
surface  row-cratering  event  are  analyzed  and  compared  to  the 
seismic  characteristics  of  single  contained  nuclear  events  to 
assess  the  quantitative  differences  between  these  two  classes  of 
events.  Finally,  the  appropriate  seismic  data  are  yield-scaled 
and  superimposed  to  simulate  the  ground  motion  associated  with  a 
PNE  evasion,  and  then  analyzed  to  determine  the  spectral  changes 
caused  by  the  clandestine  event. 

Based  on  a comparison  of  the  spectral  content  of  ground  motion 
from  the  surface  row-cratering  event,  Buggy  I,  and  the  spectral 
content  of  ground  motion  from  single-contained  events  of  equiva- 
lent yield,  and  on  an  analysis  of  the  spectral  content  of  Buggy  I 
seismograms  and  the  superposition  seismograms,  it  is  concluded 
that : 

(1)  at  all  periods  examined,  except  one,  the  spectral  ampli- 
tude content  of  single-contained  events  is,  on  the 
average,  greater  than  the  spectral  amplitude  of  Buggy  I. 
The  maximum  average  spectral  ratio  is  2.5  at  the  period 
0.5  sec(2.0Hz).  (Figure  3.5) 


(2)  in  the  case  of  scaled  seismic  data  from  the  Almendro 

event,  the  average  spectral  ratio  is  greater  than  one 
for  periods  shorter  than  0.7  seconds  (1.4Hz)  and  attains 
a maximum  of  about  2.5  (2.1)  for  a clandestine  event 

yield  of  2.5  kt  (1.375  kt)  (cf.  Figure  4.17) 

(3)  in  the  case  of  scaled  seismic  data  from  four  single- 

contained  events,  the  average  spectral  ratio  is  greater 
than  one  over  the  entire  period  range  examined.  A short 
period  maximum  occurs  at  0.15  seconds  (6.7  Hz)  and 
attains  a value  of  2.2  (1.8)  for  a clandestine  event 

yield  of  2.5  kt  (1.375  kt).  The  spectral  ratio  also  has 
a distinct  peak  at  approximately  1.0  second  (1.0  Hz) 
with  a local  maximum  value  approximately  equal  to  the 
short  period  maximum,  (cf.  Figure  4.30). 

The  above  conclusions  are  conditional,  since  they  are  based  pri- 
marily on  the  seismic  recordings  from  only  one  multiple,  surface 
row-cratering  event,  Buggy  I. 


5.2  RECOMMENDATIONS 


This  study  used  the  only  available  seismic  data  from  a surface 
row-cratering  event  to  represent  the  PNE.  In  order  to  extend  this 
study,  it  will  be  necessary  to  simulate  the  PNE  itself  by  using 
seismic  data  from  single  nuclear  device  events.  Seismic  record- 
ings are  available  from  surface  cratering  events  which  could  be 
superimposed  with  the  proper  time  delays  to  simulate  an  n-device 
surface  row-cratering  event.  With  several  simulated  PNE's  avail- 
able, a statistically  meaningful  set  of  PNE  evasion  scenarios 
could  be  examined  and  analyzed  to  define  statistical  parameters 
for  the  detection  of  a clandestine  event  in  a PNE  environment. 
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